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Subordinated processes play an important role in modeling anomalous 
diffusion-type behavior. In such models the observed constant time periods 
are described by the subordinator distribution. Therefore, on the basis of 
the observed time series, it is possible to conclude on the main properties of 
the subordinator. In this paper we analyze the anomalous diffusion models 
with three types of subordinator distribution: a— stable, tempered stable 
and gamma. We present similarities and differences between the analyzed 
processes and point at their main properties (like the behavior of moments 
or the mean squared displacement). 

PACS numbers: 05.10.Gg, 02.50-r, 02.70-c 

1. Introduction 

Brownian motion is a classical continuous-time model describing diffu- 
sion of particles in some fluid. Besides physics, it has found many real-world 
applications, like in ecology, medicine, finance and many other fields [1]. 

But in spite of many obvious advantages, the standard Brownian diffu- 
sion cannot model the real time series with apparent constant time periods 
(called also trapping events), which are often observed in datasets recorded 
within various fields. Therefore a rapid evolution of alternative models is 
observable in many areas of interest. Especially anomalous diffusion mod- 
els have found many practical applications. They were used in variety of 
physical systems, including charge carrier transport in amorphous semicon- 
ductors [SHU, transport in micelles [5], intracellular transport [6] or motion 
of mRNA molecules inside E. coli cells [?]• The constant time periods can 
be also observed in processes corresponding to stock prices or interest rates, 
so models based on the subordinated processes might be also useful in mod- 
eling financial time series, [SHlUj. 
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One of the most important issues that arises in the analysis of the sub- 
ordinated processes is the description of waiting-times that correspond to 
the periods of constant values. Finding a proper subordinator distribution 
allows to conclude on the properties of the whole process. The most popu- 
lar subordinator distribution is the inverse a— stable, see for instance 
but recent developments in this area indicate that another nonnegative in- 
finitely divisible distribution can be also used to model the observed waiting- 
times, [inillZlIS]- The family of such distributions contains, besides one- 
sided Levy stable, also Pareto, gamma, Mittag-Leffler or tempered stable. 

In this paper we analyze the subordinated Brownian motion with three 
types of the inverse subordinator distribution, namely a— stable, tempered 
stable and gamma. We show the differences between the distributions and 
present the main properties of the analyzed subordinated processes mainly 
expressed in the language of moments. Moreover, we investigate the asymp- 
totics of the mean squared displacement and show that in the gamma case 
it is linear for large t, while for small t it exhibits non-power law behavior. 



2. Subordinated Brownian motion 

We start with introducing a general definition of the considered pro- 
cesses. The subordinated Brownian motion is defined as 1131: 



Y{t) = B{S{t)), (1) 

where {B(r)}r>o is the Brownian motion and {S{t)}t>o is an inverse sub- 
ordinator of {U{T)}r>o mmS], i.e.: 

S{t) = inf{T > : U{t) > t} (2) 

for increasing Levy process {U (t)}t->o with the Laplace transform given by: 

i?e-^t^W =e-**(^). (3) 

The function ^{z) is called the Levy exponent and can be written in the 
following form 

poo 

^i>{z) = ez+ / {l-e-''')v{dx). 
Jo 

Here, 6 > is the drift parameter. If for simplicity, following [13], we assume 
6 = 0, then v{dx) is an appropriate Levy measure. Moreover, B{t) and S{t) 
are assumed to be independent. 

The probability density function (pdf) of the process {Y{t)}t>o is char- 
acterized by the generalized Fokker-Planck equation |13pi6]: 



^^^ = ^^'J>"^(x,t), u;(x,0) = 5(x), (4) 
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where 5{x) is the Dirac delta in point x and <I> - an integro-differential 
operator defined as: 

^f{t) = jj\l{t-y)f{y)dy. (5) 

The function M{t) is called the memory kernel and is defined via its Laplace 
transform, 1131: 



OO 1 

ztr - ^ 



M{z) = / e-^'M{t)dt 



3. Three cases of subordinator distribution 

The classical anomalous diffusion type model given by the subordinated 
Brownian motion ([1]) defines subordinator S{t) as an inverse a-stable pro- 
cess, see for instance |15j . It implies that the lengths of the constant time 
periods are a-stable distributed. However, in some applications also differ- 
ent distributions describing the lengths of the constant time periods might 
be useful. In this paper, besides a-stable, we consider two other cases of 
subordinator distribution, namely tempered stable TS{a, A, c) and gamma 
G{a,c). 

We start with a brief review of the main properties of the considered 
distributions. 



a-stable distribution 

Since there is no closed form for the probability density function of 
the a-stable distribution, it is usually more conveniently defined by it's 
characteristic function, given by 

i(f)-S exp-[-a"|t|" [l -i/3sign(t)tan^] if a / 1, , . 

^^ '~\exp{-a\t\[l + if3sign{t)^ln\t\]+int} if a = 1, 

where a S (0,2] is the stability parameter, /3 G is the skewness 

parameter, cr > is the scale parameter and S M is the location parameter. 
Note that if a < 1 and /3 = 1 the stable distribution becomes totally (right) 
skewed. Since the subordinator should be a non-decreasing process, in the 
following we assume a < 1, /3 = 1 and // = 0. Moreover, for simplicity we 
assume a = 1. 

Recall, that the a— stable family has two important properties. First, a 
sum of two independent a-stable random variables with the same a parame- 
ter is again a-stable distributed. Second, the tails of the stable distribution 
are governed by the power law behavior. 
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Tempered stable distribution 

The positive tempered stable random variable T with parameters q, A 
and c is defined through the Laplace transform 

E (e-"^) = e-^((^+«)'^-^"), A > 0, < a < 1, c> 0. (7) 

In the above definition A is the tempering parameter, while a and c are the 
stability and scale parameters, respectively. Observe that if A = 0, then 
the random variable T becomes simply a— stable with the scale parameter 
c^/". The probability density function (pdf) p'^^xc °^ tempered stable 
distribution with parameters a, A and c can be expressed in the following 
form: 

where o" = (c * cos 2^^^^° and p^'^^^'^'P'^) (^x) is the pdf of the a-stable distri- 
bution with the stability index a, scale parameter a, skewness /3 and shift 
[inmSllIZ] ■ Because all moments of the tempered stable distribution are 
finite, it becomes attractive in many practical applications for instance in 
finance p^fTO], biology [20] and physics like anomalous diffusion, relaxation 
phenomena [l2l[21], turbulence [22] or plasma physics [23], see also p^[25]. 

Gamma distribution 
The pdf of the gamma distribution p'-'^'^'^^ is given by 

where T{z) is the Gamma function defined as: 

/•oo 

r{z) = / t'-^e-^dt. (10) 
Jo 

It is interesting to note that for a = 1 the gamma distribution becomes 
the exponential one. Moreover, gamma distribution is infinitely divisible. 
For Xi ~ G{ai,c) we have Yll=i^i ~ Gi^'^^iOi^c) provided that Xi are 
independent. 

In Figure [3] we plot sample probability density functions, as well as, the 
tails of the considered distributions. The a parameters in stable and tem- 
pered stable distributions are equal to 0.6. The parameters of the gamma 
distribution are chosen so it's mean is equal to the mean of the tempered 
stable distribution. 



(8) 
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Fig. 1. Probability density functions (pdf) of the three considered distributions: 
S(a) - a-stable, TS(a, A,c) - tempered stable and G(a,c) - gamma. The a pa- 
rameters in the stable and tempered stable distributions are equal to 0.6. The 
parameters of the gamma distribution are chosen so it's mean is equal to the mean 
of the tempered stable distribution. Right panels display the right tails of the 
corresponding probability density functions in the double-logarithmic scale. 



4. Subordinated Brownian motion with different types of 
subordinator distributions 

In this section we examine the subordinated Brownian motion defined 
in ([1]) with three types of subordinator distribution, namely: a— stable, 
tempered stable and gamma. 

a— stable case 

The a-stable subordinator {Ua{T)}r>o is a non-decreasing Levy process 
with the Levy measure v{dx) = x~^^^'^^ Ix>odx and the following Laplace 
transform: 



Ee-^U{t) ^ ^-tz" ^ 0<a<L 



(11) 
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Therefore the function ^, that appears in Q, takes the form 

"^{z) = z". (12) 
This imphcates the form of the memory kernel, namely 

M{t) = (13) 
r(a) 

The first two moments of the subordinated Brownian motion defined in ([T]) 
in the a-stable case are given by [5]: 

<y(i)>=o, <Y\t)>= \ e, (14) 

r(a + 1) 

while the covariance function takes the form 

Tempered stable case 
The tempered-stable subordinator {Ua,x,c{T)}r>Q is a Levy process with 

tempered stable increments (i.e. with the Levy measure v{dx) = ^T+^Ix>odx) 
and the Laplace transform given by: [21] 



E (e-^^W) = e-**(^) = e-*((^+^)"-^'^), (16) 

where A>0,0<q;<1. Let us point out that in case A — >■ the operator <l> 
is proportional to the fractional Riemann-Liouville derivative, therefore ([3]) 
tends to fractional Fokker-Planck equation (see subsection a-stable case), 
[21j . The basic properties and the simulation procedure of the process 
{y(t)} defined in ([1]) in the tempered stable case one can find in |10 yi2f[H] . 

Observe that the memory kernel M{t) in the considered case can be 
calculated on the basis of the following equation 



[ e-'''M{t)dt = -f- = = — . (17) 

Jo ^(u) (A + n)" - A" ^ ^ 



As a consequence, the memory kernel M{t) takes the form: 

M{t) = e-^*t°-iS,,„((At)°), (18) 

where 

oo 

E^^p{z) = J2' 



,^„r(a^ + /3) 
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is a generalized Mittag-LefHer function, |26j . 

Since for < a < 1 and /3 < 1 + a, the generalized Mittag-Leffler 
function for z (z R and z 7^ can be expressed as (see Theorem 2.3 in |26j): 

1 r°° 

^(z) = explz^/"} + / K(a, (3, r, z)dr, 

a Jo 

where 

K(a,(3,r,z) = ^ r'^-^^/^c-^''^" "^"^""^^ " /^)) " ^ "^"(^(^ + /3 - «)) ^ 
' vra — 2rz cos(7rQ;) + ' 



(19) 

the memory kernel M{t) is given by the following formula: 

\1— Q roo 

M{t) = + e-^*r~^ / K{a, a, r, (At)°)dr. (20) 

" Jo 

Therefore we have 

lim M{t) = . (21) 

t->-oo ct 

The above limiting behavior is a simple consequence of the fact that for 
large t the generalized Mittag-Leffler function can be written as 



E^MW) = ^-^^ - t + 0(|A<|-'-). (22) 

fc=0 ^ ^ " 

Knowing the form of the memory kernel we can calculate the basic statis- 
tics of the process {Y{t)} such as moments and autocovariance function (see 
Theorem 1 in [l3]), namely 

< Y{t) >= 0, < y2(t) >= / M{u)du (23) 

Jo 

and 

/>min{t,s) 

<Y{t),Y{s)>= M{u)du. (24) 







However, in case of the tempered stable distribution such derivations require 
numerical approximations. 
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Gamma case 

The gamma subordinator {C/a,c(''")}r>o is a Levy process with indepen- 
dent gamma distributed increments, i.e. with the Levy measure v{dx) = 
c^—^Ix>odx) and the Laplace transform given by: 

Observe that in this case also the one-dimensional density p[t, x) of the 
process Ua,ciT) is given in a closed form, namely 

In this case the Levy exponent ^{z) is given by: 

^{z) = clog{l + za), (27) 
what implicates that the memory kernel M{t) can be expressed as: 

M{t) = L-^ rl , (28) 

^' \clog{l + za) J ' ^ ' 

where L~^(f(t)) is the inverse Laplace transform of the f{t) function. In 
order to find a formula for the memory kernel M(t), we use the following 
relation (being a consequence of the Proposition 1 in 



< y(t)2 >= f M{u)d 
Jo 



'u. (29) 



On the other hand, we have 



< Yity >=< S(t) >, (30) 

where {5(4)} is the inverse subordinator. Moreover, using the relation be- 
tween subordinator {?7(t)} and its inverse and the fact that for each t the 
random variable S'(t) is positive, we obtain 

/■oo noo 
< S(t) >= / P{S{t) > T)dT = / P(U{t) < t)dT. (31) 

Therefore, in case of the gamma distribution, we get 
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Table 1. Characteristics of the subordinated process {y{t)}t>o defined in ([T]) for 
the three cases of the subordinator distribution. 



Subordinator 
distribution 


Stable 
Sia) 


Tempered stable 
TS{a, A,c) 


Gamma 
G{a,c) 


Parameters 


a 


a, A, c 


a, c 






c((A + z)« - A°=) 


cIog(l + az) 


M{t) 


r(Q) 


ie-^'t«-iEc,«((At)") 


c /o aTr(T)'^^ 


< Y{tf > 


t" 
r{a+i) 


i/j*e-^"«"-i£;a,a((A-u)'^)<i« 


■JO r(cr) 


< y{t),y(s) > 


min{t,s)" 

r{a+i) 




TOO 7(cr.min(s,i)/a) 

Jo r(cT) 



where 7(,s, x) is an incomplete gamma function defined as: 

7(s,x) = / f~^e~^dt. (33) 
Jo 

Finally, from ()29p we have 

M(t = - / — - ' ' dT = / ——dr. 34 

^' cJo r(r) dt c Jo a-r(T) ^ ^ 

Again, the basic statistics of the process Y{t) can be calculated. Observe 
that from ([30|1 and ([32]) we have 

<^<"^-r^''^ (35, 

and 

< Y{t),Y{s) > = < B'^{mm{s,t) >=< S'(min(s,i) >= (36) 

7(cT, min(g,^)/a) 

^''^ 

The main characteristics of the subordinated process {Y{t)}t>o defined 
in ([1]) for the three considered cases of subordinator distribution are sum- 
marized in Table [TJ 

5. Sample paths properties of the subordinated Brownian 
motion with different types of subordinator distribution 

Sample trajectories of the process {Y{t)} defined in ([1]) are plotted in 
Figure The chosen parameters correspond to the middle panels of Figure 
[31 Observe visible differences in the character of constant time periods. 
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Fig. 2. Sample trajectories of the subordinated process {Y{t)} with the three con- 
sidered subordinator distributions: S(a) - a-stablc, TS(a, A, c) - tempered stable 
and G(a,c) - gamma. The a parameters in the stable and tempered stable distri- 
butions are equal to 0.6. The parameters of the gamma distribution are chosen so 
it's mean is equal to the mean of the tempered stable distribution. 



Now, we focus on one of the most popular characteristic of the recorded 
process trajectories in experimental analysis, namely mean squared displace- 
ment. Recall that the ensemble averaged mean squared displacement is 
defined as: 

/oo 
x'^P{x,t)dx, (38) 
-oo 

where P{x,t) is the probability of finding a particle in a infinitesimal interval 
(x, X + dx) at time t. On the other hand, the time averaged mean squared 
displacement is given by: 

c^'(t,T)= ^° rj._^ (39) 

where T is the length of the analyzed time series. 
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For a standard Brownian motion the mean squared displacement (msd) 
scales as t no matter if it is calculated as the ensemble or the time average. 
However, the behavior of the ensemble average changes under subordination 
scenario. In the a-stable case ensemble averaged msd scales as t" [27], while 
in the tempered stable case as for t — >■ and as t for t — t- oo [141121]. It 
can be shown (for a detailed derivation see Appendix) that in the gamma 
case the ensemble average scales as 

< Y\t) I * ^ °' (40) 

\ t/a as t ^ oo. 

In Figure [5] we plot the mean squared displacement calculated as the 
ensemble average over 1000 simulated trajectories in the three considered 
cases. The chosen parameters are the same as on the corresponding panels 
of Figure [3l Moreover, we fit a power law function to each of the obtained 
curves, except for small t in the gamma case. In the a-stable case the 
power law is fitted for the whole range of t, while in the tempered stable 
case separately for small and large t and in the gamma case only for large 
t. Observe that the obtained values are close to the theoretical power laws. 

Finally, we calculate the time averaged mean squared displacement. The 
obtained values calculated as the time average from a simulated trajectory 
of each of the three considered processes is plotted in Figure [5l Observe 
that in all cases the obtained msd behaves as t. 

6. Conclusions 

In this paper we have examined the anomalous diffusion models based 
on the subordinated Brownian motion with three types of subordinators 
distribution: a— stable, tempered stable and gamma. The main result is 
related to the properties of the analyzed processes. We have pointed at the 
asymptotic behavior of the mean squared displacement in three considered 
cases and showed that in gamma case for small values of the arguments we 
obtain completely different behavior (non-power) from this observed in two 
other cases. 

Appendix 

In order to show the asymptotic behavior of the ensemble mean squared 
displacement (msd) in the gamma case we use the Proposition 3 in [28] , 
namely if /i =< U{1) > is finite, then < S{t) >~ t/fi for large t, where 
{U{t)} and {^(t)} are the subordinator and its inverse defined in re- 
spectively. In the gamma case with parameters a and c, < U{1) >= ac 
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t 



Fig. 3. Mean squared displacement calculated as the ensemble average over 1000 
trajectories of the subordinated process {Y{t)} with the three considered subor- 
dinator distributions: S(q;) - a-stable, TS(q;, A, c) - tempered stable and G(a,c) 
- gamma. The a parameters in the stable and tempered stable distributions are 
equal to 0.6. The parameters of the gamma distribution are chosen so it's mean 
is equal to the mean of the tempered stable distribution. The fitted power law 
functions arc plotted with the corresponding gray lines. 



therefore when i — ?• oo we have 

< Y^{t) >=< B^{S{t)) >=< S{t) >~ — . 

ac 

In order to show the asymptotic behavior of the msd function for small t 
we use its explicit form: 



(41) 
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Fig. 4. Mean squared displacement calculated as the time average of the subordi- 
nated process {Y{t)} with the three considered subordinator distributions: S(a) - 
a-stable, TS(a, A, c) - tempered stable and G(a,c) - gamma. The a parameters 
in the stable and tempered stable distributions are equal to 0.6. The parameters 
of the gamma distribution are chosen so it's mean is equal to the mean of the 
tempered stable distribution. 



We can use the series expansion of the incomplete gamma function: 

oo 

7(5, x) = T{s)x'/^^e-'^ en(-l)xi/2"/„+.(2xi/2), 

where 



n=0 



fc=0 

and Iv{z) is the modified Bessel function defined as follows: 

/I \ " °° 



k=0 



k\T{v + k + l)' 



(42) 



(43) 



(44) 
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The function It, (2) can be for small z approximated by 

T{v + 1) 

Therefore, when t — )• 0, the function under the integral in (j4ip behaves like: 

^ ' n=0 \ ^ ^ ) 

00 



n=0 

(47) 

what gives 



7(T,t/a) (t/a)^e-*/'^ 

^X^^^ r(T + i) ■ 

Thus, when f — >■ we obtain 



(48) 



Now, note that the function can be approximated as: 

^ 00 

f77)-E«'^-' (50) 
^ ' k=\ 

for some that are independent of z and satisfy the following relation: 

n 

an = naian - 020^-1 + ^(-l)''C(^)an-fc, (51) 

fe=2 

where ^(2) is the Riemann zeta function [29]. As a consequence, we have: 

< y2(i) > V ak / (Va)"(r + l)'dT. (52) 

In order to simplify the notation denote u = t/a. We have 

p—u °° roo —u °° f-OD 

< Y\au) >~ X] "^"^^^''^^ = 5Z (53) 

^ k=i '^^ fc=i 
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Finally, let us consider the asymptotic behavior of the function /(n, k) = 
u^T^dr. Integrating by parts gives the recursive relation: 

/(n, k) = --^ - -^fiu, k - 1), /(n, 1) = --^ (l + -^"j . (54) 
log u log u log u \ log u J 

Therefore when n — )> the f{u,k) ~ what yields the asymptotic 

behavior of the < y^(an) > for small u, namely: 

< Y\au) >~ V afc = Af^, (55) 

c log u ^ log u 

where A = const. Substituting u = t/a in ([^^ we obtain: 

^ ^ clog(t/a) ^ log(t/a) ^ ^ 
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